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In a ferromagnet/heavy-metal bilayer device with strong spin Hall effect an in-plane current ex¬ 
cites magnetic dynamics through spin torque. We analyze bilayers with perpendicular magnetization 
and calculate three-dimensional phase diagrams describing switching by external magnetic field at 
a fixed current. We then concentrate on the case of a field applied in the plane formed by the film 
normal and the current direction. Here we analytically study the evolution of both the conven¬ 
tional “up”/“down” magnetic equilibria and the additional equilibria created by the spin torque. 
Expressions for the stability regions of all equilibria are derived, and the nature of switching at each 
critical boundary is discussed. The qualitative picture obtained this way predicts complex hystere¬ 
sis patterns that should occur in bilayers. By analyzing the phase portraits of the system we show 
that when the spin torque induced equilibrium exists, switching between “up” and “down” states 
proceeds through it as an intermediate state. Using numeric simulations we analyze the switching 
time and compare it to that of a conventional spin torque device with collinear magnetizations of 
the polarizer and the free layer. 


I. INTRODUCTION 

Recently a number of investigations focused on bilayer 
structures consisting of a ferromagnetic (F) layer and a 
nonmagnetic (N) layer with strong spin-orbit (SO) inter¬ 
action made of heavy-metals such as Pt, Ta or W.^^^^^It 
was theoretically predicted and experimentally observed 
that when an in-plane electric current is being applied, 
the itinerant electrons inside the nonmagnetic layers be¬ 
come spin polarized due to the strong spin-orbit coupling 
and exert a spin torque on the ferromagnetic layers. This 
additional torque contributes to the magnetization dy¬ 
namics described by the Landau-Lifshitz-Gilbert (LEG) 
equation. Up to now two different models have been pro¬ 
posed to account for the effects. One of thenJ^ treats 
the bilayer structure as a two-dimensional system with 
strong interfacial Rashba spin-orbital coupling due to the 
structural inversion symmetry breaking in the direction 
normal to the interface.^ This model leads to a field-like 
torque directed along mx (jg where m = M/M^ 

is the magnetization unit vector (Mg represents the con¬ 
stant absolute value of the magnetization M), jg is the 
in-plane electric current density, and the interface nor¬ 
mal v ector is pointing along the z direction. The other 
modePEEHU ig based on the interfacial diffusion of the 
pure spin current that originates in the heavy-metal lay¬ 
ers due to the bulk spin Hall effect (SHE]^J^ and leads 
to spin transfer torque (STT) dynamic j^^ l ^^ * in the mag¬ 
netic layers. In the SHE model the torque is directed 
along ih X [ih X (jg X z)].l^This type of torque is called a 
Slonczewski, or damping-like, or adiabatic torque in the 
literature. 

Several experiments showed that an in-plane electric 
current flowing through the structure is able to switch the 
magnetization of the ferromagnetic layerIn those 
experiments the F layers were magnetized perpendicular 
to the film plane. It is believed that the observed mag¬ 


netic reversal can only be induced by the damping-like 
torque. The reasons for this are (a) experimentally mea¬ 
sured switching phase diagr ams are in accord with the 
macrospin model calculationJ^I^, and (b) due to its sym¬ 
metry, a field-like torque, if there is any, does not favor ei¬ 
ther “up” or “down” state of the perpendicular magneti¬ 
zation, and therefore should not contribute to switching. 
These arguments seem to favor the SHE-based model, 
however, subsequent calculation^^^M^ suggested that the 
model based on Rashba coupling generates both field¬ 
like and damping-like torques, and thus is also capable 
of describing the switching. 

Despite the fact that the underlying torque mecha¬ 
nism is still not fully understood, a thorough study of the 
switching behavior based on the existing experimental re¬ 
sults is of importance for analysis and prediction. In this 
paper we study magnetic switching induced by an exter¬ 
nally applied magnetic field H at a fixed in-plane electric 
current. It is assumed that the magnetic anisotropy en¬ 
ergy of the F layer is uniaxial 

F^(m) = —iF(m • z)^. (I) 


z 



FIG. I. Schematic diagram of the bilayer spin Hall effect 
device. Electric current j flows along the y axis. 
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Magnetization is switched between “up” and “down” 
at critical fields He that form a surface in the three- 
dimensional i7-space. Without electric current and for 
the magnetic anisotropy given by 0 this surface is an 
axially symmetric figure of revolution with a cross-section 
given by the astroid curveP^The presence of the in-plane 
current je breaks the axial symmetry of the critical sur¬ 
face He- To find its shape we use the method of Refs. [29] 
and EOl A related problem of finding the critical field 
magnitude as a function of current Hc{je) for a given 
field direction was considered numerically in a number 
of publicatoins.l^ In our treatment we consider only the 
damping-like torque. Since the field-like torque can be 
compensated by an external field, its presence simply 
shifts the critical surface and we do not include it in 
the calculations. First, we find a three-dimensional crit¬ 
ical surface. Second, we study a particular section of 
this surface with H confined to the plane formed by the 
directions of current je and magnetic easy axis. Here 
we derive analytic formulas for the switching bound¬ 
aries He. We also observe that, in accord with numeric 
invest igations,!^ in the presence of a sufficiently large cur¬ 
rent an extra equilibrium direction appears in addition to 
the “up” and “down” equilibria, and get analytic expres¬ 
sion for its location. By analyzing the phase portraits of 
the system we show that if additional equilibrium exists, 
switching form “up” to “ down” state proceeds through 
this intermediate state. Using numeric simulations we 
analyze the switching time and compare it to that of a 
perpendicular polarizer device. 


II. GENERAL DESCRIPTION OF THE 
THEORETICAL APPROACH 


Magnetization dynamics of the ferromagnetic layer in 
the macrospin approximation is governed by the LLG 
equation: 


dfVI ^ , a 

-dt- =“•«>+ jf. 



( 2 ) 


where a is the Gilbert damping factor, 7 is the gyro- 
magnetic ratio, and Heff is the total effective field. The 
standard LLG equation can be transformed into 


—— = — ih X heff — am x ih x heff, 
at 


( 3 ) 


where the field is rescaled as heff = Heff/i^/c using the 
characteristic anisotropy field Hj^ = 2K/fioMs^ and the 
time is rescaled as t ^ (1 a‘^). Hereafter, 

all the field-related terms that are written in lowercase 
letters are di mensio nless (normalized by Hk). 

The methoeP^^ can be summarized as follows. A sta¬ 
tionary solution iho of Eq. ^ satisfies the equilibrium 
condition ih x heff = 0, which indicates that the mag¬ 
netization at equilibrium should be parallel to the total 
effective field, i.e., heff = Aiho with arbitrary A. Total 


effective field is given by heff = h — Ve + hgp, where h 
is the external field, 5 = E/{fioHkMs) is the rescaled 
anisotropy energy, and hgp is the spin torque effective 
field 


hsp = aj[m X (e^ X 5)], (4) 

where ej is a unit vector in the electric current direction 
and aj is a spin torque strength parameter, proportional 
to the electric current density. Equation heff = Aiho can 
be solved for the external field as h = Amo + \/e{mo) — 
hgp (mo). The meaning of this formula is that for any 
given magnetization direction there is a whole line of ex¬ 
ternal fields, parameterized by A, for all of which this 
direction is an equilibrium (stable or unstable). In spheri¬ 
cal coordinates with three orthogonal unit vectors defined 
as ih = sin6>cos0x + sin6>sin0y + cos6>z, 0 = dm/dO^ 
and (j) = {1/ sin^)9m/90, we get 

h = Aiho + {dee - Oq + , (5) 

where de stands for d/dO and the superscript 0 indicates 
the 0 component of a vector field (e.g., = heff • 6>), 

etc. Equation ^ maps the 3D space {A,6>o,0o} to the 
3D field space {TTa,, hy^hz}. 

Next, stability of equilibrium states is analyzed. This 
is done by expanding Eq. 0 in small deviations ih = 
iho + (5m up to the linear terms, such an expansion pro¬ 
duces two coupled linear differential equations 

with matrix A (^05 0o) given by 

^ ^ Ideiah^g + 7ff) ^d^iah^g + hfg)'] 
[deiahfg - h^g) ^d^{ah*g - hlg)\ 

Stationary solutions can be classified as stable or unsta¬ 
ble using the eigenvalues of the matrix A, which are 
uniquely determined by its determinant, detA, and trace, 
trApSl An equilibrium is stable when both eigenvalues ii± 
are either complex numbers with negative real parts, or 
negative real numbers, which leads to the stability crite¬ 
ria 


detA > 0 and trA < 0 , ( 8 ) 


that are applied to select the subset of the line h(A, 6 >o, 0o) 
for which (po) equilibrium is a stable, i.e., to find the 
values of A for which conditions ^ are satisfied. By eval¬ 
uating expression Q at external field specified by Eq. ^ 
one obtains A(A). We find that for the arbitrary form of 
spin torque and arbitrary anisotropy energy, trA (A) is 
linear in A with a negative coefficient, and detA (A) is a 
quadratic function of A with a positive quadratic coeffi¬ 
cient. To simplify the expressions, we introduce a vector 
field f = — + hgp and its gradient 


Vf = 


def 

d^f 


( 9 ) 
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(see Appendix [A| for the explicit expressions). The roots 
of equations trA(A) = 0 and detA(A) = 0 can then be 
respectively calculated as 


At(6^o, 0o) = 2 


def + 84 , f ^ {def - d^f) 


^±{(^ 0 , 0 o) — 


def + d^f4> 


( 10 ) 

( 11 ) 


± 




In terms of the critical values At and A±, the stability 
criteria become 


[a > Max(AT,AT) if At ^ A_, 

I At < A < A_ or A > At if At < A_. 


( 12 ) 


When At are complex, detA is always positive and cri¬ 
teria ( 12 ) can be further simplified as A > At- The full 
classification is given in Appendix |B| 

Substituting functions A^, A_ or At for A in Eq. (§ 
one generates three surfaces in the field space, which we 
denote as 5'+, S- and 5 't respectively. Their physical 
meaning is that at least one equilibrium changes its sta¬ 
bility when external field crosses such a surface. The 
equilibrium is either locally destabilized when the 5 't 
surface is crossed, or merges with a saddle point when 
the S± surfaces are crossed.!^ The entire critical surface 
S is constructed from the parts of 5'+, S- and 5 't as 
explained in Ref. [30l 


III. 3D PHASE DIAGRAM 




<P 


FIG. 2. (color online) The critical A surfaces (upper panel) 
and the cross-section at 0 — 7r/2 (lower panel). The A+, 
A- and the At surfaces are plotted in blue, green and red, 
respectively. The A- surface is always below A+. The At 
surface in this case is separated into two parts, one above 
A+ and the other below A_. In the lower panel, stable and 
unstable regions are denoted by white and gray, separated by 
the active parts of the three lines. The damping parameter 
and the spin torque strength are set to be o = 0.1 and oj = 
0.1. 


In this section we construct the three-dimensional crit¬ 
ical surface using the method of Sec. [Ilj The dimen¬ 
sionless perpendicular anisotropy energy has the form 
5 = — cos^ 0/2. We set the in-plane current to be along 
the +y direction, and the current induced field is then 
hgp = Of jin X X with aj given by aj = 6 >sHie/io 5 where 
^SH is the spin Hall angle, jo = 2 eM^dYH]^/h is the char¬ 
acteristic current density, and dp is the thickness of the F 
layer. For brevity, we drop index “0” for the equilibrium 
direction. The critical values of A calculated according 
to ( 10 ) and specialize to 


Xt{0, </)) = -cos^ 0 — — sind cos 0, (13) 

2 a 


A±(d,j)) = 


sin^ 9 


' 0 ±i 


I sin^ 0 


— a‘j cos^ ( 


(14) 

For aj = 0 one finds At to be in the midpoint of the 
interval (A_, A^) for any direction (d,0). Then, accord¬ 
ing to the criteria ( 12 ), 5'+ constitutes the entire critical 


surface. For non-zero current At may leave the inter¬ 
val (A_, At) for some values of (d,0). In those cases 
switching boundaries should be selected for every direc¬ 
tion individually. Fig. shows the three critical values as 
functions of equilibrium angles at nonzero current with 
aj = 0.1 and a = 0 . 1 . 

For practical calculations we decompose Eq. (§ into 
Cartesian coordinates 

= (A + cos^ d) sind cos 0 , (15a) 

hy = {X-\- cos^ 0) sin d sin 0 — aj cos d, (15b) 

hz = {X — sin^ 0 ) cos 0 + aj sin 0 sin (j). (15c) 

Eor arbitrary A’s these equations represent the mapping 
of the (X^O^cj)) space to the {hx^hy^hz) space for the 
case of uniaxial anisotropy and chosen electric current 
direction. If functions Xt{0,(I)) @ or X±{0,(f>) @ are 
substituted for A, one obtains parametric expressions for 
the critical surfaces St and S± with parameters (d, 0) 
running through all possible values, 0 ^ ^ tt and 

0 ^ ^ 27r. The 3D phase diagram in field space with 
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FIG. 3. (color online) The critical surface in three- 
dimensional field space. The same color convention and pa¬ 
rameters as in Fig. are adopted. 


the same parameters is shown in Fig. 


IV. PHASE DIAGRAM IN THE Y-Z PLANE 


The 3D phase diagram is quite difficult to use due to 
the complicated shape of the critical surface S. The re¬ 
lated experiments are often perfor med with the field be¬ 
ing confined within the y — z plane.^^^^ Here we study in 
detail a section of S corresponding to the external field 
h confined to such a plane. This section is a line S in the 
2D space {hy, hz). 

A field in the y — z pl ane satisfies a constraint hx = 0. 
According to Eq. (15a) this implies a relationship be¬ 


tween 0^ 0, and A. On the one hand, this relationship al¬ 
lows one to express the equilibrium angles ( 6 >, 0 ) as func¬ 
tions of {hy^hz) and study how the equilibria evolve as a 
function of applied field. On the other hand, Eq. (15a) 


can be used to find the section S. While the surface 
S is given by a parametric formula w ith independently 
varying 0 and 0 as explained in Sec. the line S is 
found from the same formula with 0 and 0 be ing related 
to each other by a requirement that Eq. (15a) holds with 
A = At,+,-(6^, 0). 


A. Evolution of equilibrium states 


Equation (15a) has three types of solutions: (I) (j) = 
:^7rl2, (II) A = — cos^ 6 >, and (III) = 0 , 7 r. Since the 
value of 0 at 6 > = 0 , 7 r is immaterial, type III can be 
considered as a special case of type I. Thus we focus on 
the first two cases. Eor definiteness, assume aj >0. 



FIG. 4. Equilibrium states shown as points on the unit sphere 
for small aj > 0 in the absence of external field. Arrows show 
the directions of equilibria displacements as an external field 
h is increased at constant aj. Points OMi, OM 2 , and S 2 
eventually merge into X. 


Solutions of type I have sin(/) = ± 1 . They are located 
on the meridian of the unit sphere lying in the y — z plane. 
We w ill ca ll the m on- meridian states. Eliminating A from 
Eqs. (15b) and (15c) one finds a system of equations for 
their polar angles 


<p = X'K j2 , 

hy COS 0 hz sin 0 = X sin 0 cos 0 — aj. (16) 


Depending on hy, hz and aj there can be four, two or 
zero equilibrium states of type I. 

Sol ution s of type II have A = —cos‘^0. Equations 
(15b) and (15c) read 


hy = —aj cosO, 

hx = — cos 0 + aj sin 0 sin (j). 


Solving them one finds 


cos 6 > = 



sincj) = 


h>y j aj 



(17) 


with associated requirements \hy\ < aj and \{hz — 
hy/aj')/ — hy\ < 1 that define their domain of ex¬ 
istence. Having 0 7 ^ ffi 7 r/ 2 , type H solutions are away 
from the y — z plane and will be called off-meridian states. 
They exist as a pair with the same polar angle 0 and az¬ 
imuthal angles (j) and tt — (j). 

Equilibrium states can be visualized as points on the 
unit sphere that change their positions when the experi¬ 
mental parameters h and aj are varied (Eig. |^. In the 
absence of current and external field the uniaxial mag¬ 
net exhibits two stable equilibria 0 = 0 , tt: the “up” and 
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“down” states. Due to the axial symmetry of the sys¬ 
tem the entire equator of the unit sphere forms a circle 
of unstable equilibrium state. 

At non-zero current spin torque breaks the axial sym¬ 
metry of the problem even in the absence of magnetic 
field. For 7 ^ 0, h = 0 the continuous set of unsta¬ 
ble equilibria along the sphere’s equator is reduced to 
four isolated equilibrium points. Two of them are off- 
meridian states (OMi ^2 in Fig.[^ for which Eq. ( [I7| ) gives 
0 = 7 r /2 and 0 = 0 , 7 r, i.e., the directions. The other 
two are on-meridian states Si ^2 given by Eq. ([T^. It 
will be shown below that they are saddle points. Eor 
small aj the system has six equilibria: slightly displaced 
“up” (Mup) and “down” {Mdown) on-meridian states, on- 
meridian states Si ^2 that are slightly displaced from the 
equator of the sphere, and the states OMi ^2 (Fig.®. 
As aj is increased, the states Si and S 2 approach the 
“up” and “down” states respectively. At a critical cur¬ 
rent they merge pairwise and disappear. 

The following useful rules apply to the on-meridian 
equilibria described by Eq. ( [l^ : (1), Increasing current 
shifts points clockwise and points Si ^2 coun¬ 

terclockwise along the meridian; (2), Increasing magnetic 
field shifts points Mupidown along the meridian towards 
the field directions and points Si ^2 away from the field 
direction (see Eig. [^. 

Consider now the situation with a small fixed current 
and a variable external field. Eor the discussion we will 
assume a fixed direction of h between -^y and Sz di¬ 
rections (see Eig. [^. Equations 0 show that as the 
field magnitude h is increased, the off-meridian states 
approach the meridian and reach it at a critical field 
magnitude. Since the two off-meridian states are mir¬ 
ror symmetric with respect to the y — z plane, they reach 
the meridian simultaneously and merge. Moreover, using 
Eq. ( [TtI ) one can show that the merging point also sat¬ 
isfies Eq. ( p!^ , so actually a merging of two off-meridian 
and one on-meridian equilibrium takes place. This tri¬ 
equilibrium merging is not accidental. As discussed in 
Ref. |34l merging of equilibria has to conserve the wind¬ 
ing number and it would be impossible for the two off- 
meridian equilibria with equal winding numbers to merge 
without the participation of a third equilibrium with the 
opposite winding number. 

As h is increased further, equilibria S 2 ^ OMi, and OM 2 
merge into one equilibrium X that remains on-meridian. 
Analysis in the next section shows that X is an unstable 
focus, analogous to the maximum energy point equilib¬ 
rium of a conventional (no spin torque) uniaxial magnet 
subjected the external field. In general, above the critical 
field the evolution of the four on-meridian equilibria M^p^ 
^downi Si^ and X is qualitatively similar to that found 
at aj = 0. We may conclude that our system has two 
qualitative regimes: one at low magnetic field where spin 
torque dominates, and another one at high field where 
magnetic torque dominates. The spin torque dominated 
regime is characterized by the presence of two OM equi¬ 
libria produced by current. In the field dominated regime 


the current-induced equilibria are gone. 

These results are quite natural. The SHE system is 
equivalent to a conventional spin-transfer device with 
spin polarizer p directed along -\-x. Spin torque attracts 
the magnetization to p and repels from — p. At very 
large currents spin torque dominates all other torques, so, 
only two equilibrium points — one close p and another 
close to — p should remain. In our system we find that it 
is enough for the spin torque to dominate the magnetic 
field torque in order to produce these equilibria. This 
happens because for the easy axis energy 0 and p = X 
the anisotropy torque is equal to zero at m = ±p. 


B. Stability of equilibria analysis and switching 
phase diagram 

In this section we are going to find the critical line S 
of equilibrium destabilization. It will be composed from 
parts produced by type I and type II solutions. 


1 . On-meridian equilibria 


Equations (11^ and (14) show that for the on-meridian 


states At is the midpoint of A± interval for any current 
value. Therefore only A^ is needed to calculate the criti¬ 
cal surface . By substi tutin g A = (j)) and sin0 = ±1 

into Eqs. and ( |15c[ ), we get an exact parametric 

form of aS'm, the line of on-meridian equilibria destabi¬ 
lization. It turns out to be the same as the one found 
befor^^ using an approximate method. 


hy = ± sim 0 — aj cos 0 , 
hz = — cos^ 0 ± aj sin 0 . 


(18) 


By evaluating det A and trA for each on-meridian equilib¬ 
rium it is possible to conclude that the “up” and “down” 
equilibria are stable foci, while the Si ^2 equilibria are 
unstable saddle points. The Sm curve for various spin 
torque strengths are shown in Eig. 

When magnetic field h crosses Sm^ one of the sta¬ 
ble equilibria M^pjdown merges with one of the sad¬ 
dles Si ^2 and disappears. In fact, Sm represents the 
Stoner-Wohlfarth (SW) astroid boundary modified by 
the current!^ The original astroid shape is squeezed 
along one of the diagonal directions. The segments of Sm 
connecting the corner points of the astroid become un¬ 
equal: two of them grow with increasing and the other 
two shrink. The details of merging depend on whether 
the long or the short segment is crossed by the field. On 
the short segment the sign of the rriz component of the 
disappearing equilibrium is always opposite to the sign of 
the field component hz^ This property was satisfied ev¬ 
erywhere on the conventional SW astroid boundary, and 
we denote the short segment of Sm as Smc with index 
“c” meaning “conventional”. On the long segment the 
sign of VTiz is not determined by the sign of hz^ Indeed, 
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FIG. 5. (color online) Lines Sm (black) and the critical field 
boundary where OMi, OM 2 and aS '2 merge into X (red/grey). 
The spin torque strength aj is set to be a. 0.0, b. 0.3, c. 0.5, 
and d. 0.7. 

points Mup and Si merging on this segment have rriz > 0, 
and at the same time it can be crossed by a field with 
> 0, if h is directed at a small enough angle to the 
y axis. We denote the long segment as Smu with index 
“u” meaning “unconventional”. 


- 0.5 0 0.5 



FIG. 6. (color online) The Bom boundaries for various damp¬ 
ing factors with spin torque strength set to aj — 0.5. The 
damping parameters a/aj are set to be 0.8, 1.5, 1.9, and 2.1 
(going from the outermost to the innermost curves). At low 
damping regime. Sqm is approaching the tri-equilibria merg¬ 
ing boundary which is shown as a red (grey) dashed curve in 
the figure. 


2. Off-meridian equlihria 


Next, we analyze the stability of the off-meridian equi¬ 
libria. According to Eq. (14) A ^ A_ is automatically 


satisfied when A± are real. Thus, according to criteria 
( [T^ , only the St critical line is relevant whether A± are 
real or complex, and the stability condition is given by 


A — At = ( ^ cos (f) -sin 0 ) sin ^ > 0. 

a 2 


(19) 


This inequality can be satisfied only for cos0 > 0, 
which means that the OM 2 equilibrium characterized by 
7r/2 < (f> < 37r/2 is always unstable. The off-meridian 
equilibrium with —it/ 2 < <f < 7 t/2 can be stable. The 
critical line A = At gives a destabilization boundary Sqm 
for t his eq uilibr ium. Its analytic form is obtained from 
Eqs. (|15b|) and (|15c|) as 


a I 


a,. 2 - fty 


ll- 


4a 


d n 


1 - (20) 


The SoM curve for various damping parameters at a fixed 
spin torque strength is shown in Eig. Erom det A and 
trA analysis it is possible to extract more detailed infor¬ 
mation about the nature of the OMi ^2 equilibria. We find 
that OM 2 is always an unstable node (two real positive 
eigenvalues), and OMi is a stable focus (complex conju¬ 
gate eigenvalues with negative real parts) inside its do¬ 
main of stability bounded by the Sqm line (see Appendix 


[C] for a complete analysis). As the field increases and 
moves out of this domain, OMi is destabilized but not 
destroyed. It continues to exist, first as an unstable fo¬ 
cus, and then as an unstable node, until it finally merges 
with the points OM 2 and S 2 ^ as discussed in Sec. iVAl 
More details of the OMi state evolution are provided in 
Appendix 0 


C. Current-field diagrams 


Switching phase diagrams Eig. and Eig. describe 
experiments performed at fixed current with magnetic 
field of a fixed direction increased until switching hap¬ 
pens at a critical value he- In a more complicated ex¬ 
periment one can measure how the he threshold depends 
on the current magnitude. Such experiments were in¬ 
deed recently performed by Yu et al^ and investigated 
numerically.!^ The hc{dj) dependencies were measured 
for different field directions and, quite surprisingly, it was 
found that for fields making a finite angle with the y axis 
multiple switchings may occur. This fact finds a natural 
explanation in the framework of our theory. In terms of 
Eig.j^the critical fields are determined by intersections of 
a straight line representing fields of given direction with 
the lines Sm and Sqm- If the direction of the field is 
defined by the angle Oh with the z-axis, the former inter¬ 
section point hcM{oLj) can be found by solving Eqs. (18) 
with hy = /icsin^/i, hz = he cos Oh- Eor the latter inter- 
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FIG. 7. (color online) Tilted current-field phase diagram at 
a = 0. The black and red (gray) lines correspond to hcM and 
hcOM^ respectively. The solid lines represent the diagrams 
with a tilting angle of the field set to be 15 degree with respect 
to y direction. The dashed lines represent the diagram with 
h\\y. Opposite tilting happens when the angle is negative. 


section point hcOM{(^j) Eqs. ( [^ should be used. The 
results are shown in Fig. One can see that hcM{oLj) 
exhibits a sharp peak located at aj = 0 for Oy = 7r/2. 
As the field is tilted away from the ^-axis, the position 
of the peak moves and its initially symmetric shape de¬ 
forms. Eventually the deformation grows so big that the 
function hcM{(^j) becomes multi-valued, in accord with 
experimental findings. Comparing the current-field dia¬ 
gram with the experimental diagram (Fig. 3 of Ref. [2^ 
one can see a good qualitative correspondence. 

Here we show how the shape of hcM{(^j) can be un¬ 
derstood from the evolution of the modified astroid Sm- 
As the current is increased with aj > 0, the astroid is 
squeezed in the (1,-1) direction and expanded in the 
(1,1) direction. The Smu lines approach the origin in 
the /i-space and the Smc lines move away from it. When 
the field is directed along the y axis, its line intersects 
the Smu boundary. Since this boundary moves towards 
the origin with increasing the function hcM{(^j) is de¬ 
creasing. However, when the field is directed at an angle 
to the ^-axis, its line may initially cross the Smc bound¬ 
ary. Since Smc nioves away from the origin, the function 
hcM{oLj) would increase. At a threshold value of current 
the field line goes exactly through the corner point be¬ 
tween the Smc and Smu- At this point hcM{oLj) exhibits 
a cusp. For currents above the threshold, the field line 
crosses Smu and, just like in the h||^ case, hcM{oLj) be¬ 
comes a decreasing function. For some angles Oy there 
may be situations when the field line crosses both Smc 
and Smu lines. This is when hcM{(^j) becomes multi¬ 
valued and complicated hysteresis patterns are realized. 

The form of the other critical field line, hcouio^j) (red 
curves in Fig. [^, is related to the evolution of the Sqm 
line. Since this line moves away from the origin in all 
directions, hcOM{o^j) turns out to be an increasing func¬ 
tion. 


D. Discussion of the phase diagram 


The Sm and Sqm lines together give the complete 
switching phase diagram in the y — z field plane. For 
small values of aj the critical line Sm is qualitatively 
equivalent to the conventional Stoner-Wohlfarth astroid, 
and the equilibrium merging process is similar: There are 
four on-meridian equilibria for h inside the astroid, and 
as the field crosses its boundary two of them merge and 
disappear. Above the critical current ^ = 1/2, the Sm 
critical line becomes self-crossing (Fig.|^. At the critical 
current the Smu lines touch each other at h = 0 so the 


threshold can be found from Eq. (18) with hy = hz = 0. 


Inside the region of self-crossing there are no on-meridian 
equilibria, as already observed in Ref. [31] and the supple¬ 
mental material of Ref. Ul However, the Poincare-Hopf 
theorem is not violated due to the presence of the off- 
meridian equilibria. 

In the absence of current the system in constant exter¬ 
nal field h resides in one of the two stable M equilibria. 
As aj is increased, the oval-shaped region of stability of 
the OMi state grows, and the area inside the modified 
astroid Sm shrinks. Moreover, the self-crossing region of 
Am, where no on-meridian equilibria exist, also grows. 
Thus both Mup and Myown states eventually become un¬ 
stable at some critical current a^ and m switches to the 
OMi state. What happens if the current is subsequently 
decreased? The answer to this question can be read from 
the hc{aj) dependence shown in Eig. At a given h 

le do^ 


the off-meridian state remains stable down to the cur¬ 
rent a^^ obtained from the equation h = hcOM{otj). If 
a^^ < a^ ^ one would observe hysteretic behavior of the 
system in the current interval between a^^ and a^ ^ At 
the higher end of this interval the system switches from 
an M state to the OMi state. At the lower end it switches 
back to an M state. As seen from Eig. the length of 
the bistable interval becomes larger for smaller h. At 
h = 0 one finds using Eq. (20) that OMi is stable for 


aj > a/2. The higher end of the interval was already dis¬ 
cussed — it corresponds to the first self-crossing of Am, 
i.e., to aj = 1/2. Eor typical values of Gilbert damping 
Of ^ 0.01 the hysteresis range is very large. It requires an 
initial pulse of current of the order aj ^ 1 to get to the 
OM state, but after that the current can be reduced to 
aj ^ (a, and the OM state can be comfortably studied 
at low currents. Experiments with SHE deviceare 
already performed in the regime aj ^ 1 so the discussed 
hysteresis should be observable. 


When magnetic field is set inside the domain of exis¬ 
tence of OM states but outside of their domain of sta¬ 
bility, the system has two unstable OM equilibria. It 
is possible to arrange parameters so that there no M 
equilibria either (this happens in the high damping, high 
current regime). In this case the system has no choice 
but to follow some precession cycle, the analysis of which 
is beyond the scope of the present paper. 








V. DYNAMIC PROPERTIES 


B. Two-stage switching through the ofF-meridian 
equilibrium 


In this section we discuss what happens after the sta¬ 
bility boundaries are crossed and equilibria are destabi¬ 
lized. 


A. Switching to the ofF-meridian state 


Existence of a stable OM state within the area given 
by Eqs. (20) raises a question: When an M state is desta¬ 
bilized at the Sm boundary, will the system switch to the 
other M state, or to a stable OM state? To answer this 
question we plot the flow diagrams (phase portraits) of 
the system. The results of calculations are presented in 
the form of qualitative sketches that emphasize the struc¬ 
ture of the flow (Eig. [^. In the field-dominated regime 
the flow is qualitatively similar to that in the absence of 
the current. There are two basins of attraction of stable 
points Mup and Mdown (white and dark areas in Eig.l^). 
The separatrix between the two basins winds around the 
unstable focus X making infinite number of turns. As 
a result, near X the basins are finely intermixed and a 
small change in initial conditions may change the equilib¬ 
rium where the system ends up. When modified astroid 
boundary is crossed, one of M points is destroyed. A 
system initially residing in this point will switch to the 
other M point. 


In the current-dominated regime there are three basins 
of attraction. The one of the OMi point (dark area in 
Eig. 1^) separates those of Mup and Mdown (white ar¬ 
eas). The white areas touch at the point of unstable 
equilibrium OM 2 . The important difference from the 
field-dominated regime is that OM 2 is an unstable node, 
and not a focus. Thus, there is no winding of the sep¬ 
arating line around it and, consequently, no area of fine 
intermixing. Eig.|^ shows what happens when Mup and 
Si points merge at the modified astroid boundary. The 
phase portrait in the upper part of the unit sphere qual¬ 
itatively changes: The basins of attraction of Mup and 
OMi merge, forming one bigger basin of attraction of 
OMi. This transformation of the phase portrait does 
not affect the qualitative picture in the lower part of the 
unit sphere and the boundary between the basin of OMi 
and Mdown- The end result is that a system initially re¬ 
siding in Mup will switch to the OMi state with certainty. 
The latter statement, of course, only applies to the case of 
slow, quasi-static change of parameters, in which case one 
can be sure that m follows the stable point with great 
accuracy. If parameters are changed at a finite speed, 
there will be a lag between m and the equilibrium point, 
and a more careful investigation would be required. 


Magnetization reversal is one of the most important 
processes in magnetism that is linked to the magnetic 
data storage process, such as in hard disk drives. Switch¬ 
ing speed and reliability are two crucial factors to the de¬ 
sign of such systems. In conventional spin-transfer torque 
switching spin polarizer is directed along the easy axis of 
the free layer. Then the magnetization moves towards the 
new equilibrium al ong a spi ral trajectory in a reliable but 
fairly slow manner faster reversals, which 
are often called precessional switchings, have been de¬ 
signed. Some have magnetic field applied orthogonally 
to the easy axis. Others use spin polarizer perpendicu¬ 
lar to the easy plane of the free layer (“magnetic fan” 
geometry ) both cases the reversal processes is 

fast but requires current or field pulse length to be per¬ 
fectly adjusted. This is experimentally hard to control 
and increases the error rate. 


It was numerically found in Refs. |32] and [35] that 
switching from an M state to the OM state is fast and 
does not exhibit precession. A recent micromagnetic sim¬ 
ulation paper has reported the similar switching behav¬ 
iors, which indicates the validity of the OOP switching 
in the microspin regime .33. Eigure shows the process 
of switching from Mup to OMi state. It is seen that the 
switching time is of the order of ferromagnetic resonance 
period T {T = 2ti in dimensionless unit used in Eig. [^. 
We will consider one possible scenario of switching from 
Mup to Mdown equilibrium with an intermediate stop in 
the OMi state. Consider a system that is initially in the 
Mup state. Magnetic field h is set in the negative 2: direc¬ 
tion during the whole switching procedure with h < 1 so 
that Mup = is stable. Eirst, we apply a short pulse of 
strong current aj ^ 1. The rise and fall times of the pulse 
are assumed to be negligible. During the pulse time Mup 
does not exist and m switches to OMi . The pulse length 
tp is selected to be large enough for the switching pro¬ 
cess bo be accomplished. Importantly, this requirement 
sets only a low bound for tp — there will be no harm 
in keeping the current on for a longer time. According 


to Eqs. (20) for hy = the state OMi has 6> = 7r/2 and 
sits on the equator of the unit sphere. After the end of 
the pulse the current is switched off and the second stage 
of switching begins. Now the states Mup to Mdown are 
stable again and m should go to one of them. With field 
pointing down and aj = 0, the boundary between the 
basins of attraction of Mup and Mdown is a parallel cir¬ 
cle, located above the equator of the unit sphere. Thus 
the second stage starts with m residing in the basin of 
attraction of Mdown^ to which m eventually relaxes in a 
precessional manner. The whole process is characterized 
by a fast first stage with strong current and a slow second 
stage, during which the systems is not driven externally. 
While the total switching time is of the same order of 
magnitude as in the conventional switching, the “active” 
stage requires much shorter time, comparable to that of 
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FIG. 8. (color online) Sketches of the phase portraits on the stereographic plane. The dashed circle in each subfigure is the 
projection of the equator. The north and south poles are projected to the circle center and to inhnity, respectively 



FIG. 9. (color online) Evolution of the three components 
of the magnetization in a two stage switching process. The 
parameters are: hy = 0.15, aj = 0.45, and a — 0.05. The 
held-current pulse is turned on at t = 0 and off at t = 20. 


precessional switching, making the procedure potentially 
useful for special applications. An important drawback 
of this switching scenario is that for a given direction of 
h it can be performed only in one direction, e.g., in the 
discussion above from Mup to Mdown- To switch back 
one would have to reverse the direction of h. 

It is interesting to to note that the SHE device switch¬ 
ing between M^p and Motown a two-stage manner de¬ 
scribed above can be alternatively viewed as a realization 
of a controlled-NOT gate with hz being the control pa¬ 
rameter. 

Finally, we want to compare the duration of the fast 
stage of SHE switching with the switching time of a con¬ 
ventional spin-torque device, where the magnetic polar¬ 
izer and the external field are both pointing along the 
easy axis of the free layer. Assuming the conventional 
spin torque to be described by a constant spin-transfer 


efficiency factor g{0) = g, we get = aj[m x z] 
with aj = gj/jo for its effective field. In this fully 
axially-symmetric case the switching time can be com¬ 
puted analytical!}®! as 

^ 2a{l rriMlh'') ’ 

with h' = h -\- aj/a and rrizo being the initial value of 
the magnetization component along the easy axis. In 
the small-damping (a <C 1), large-current {otj/a ^ ^ 

1), regime this simplifies to ~ — ln(l — rriz^^j{^aj\ 
Conventional switching requires some initial deviation of 
m from equilibrium. This deviation is usually thought 
to come from thermal fluctuations and can be evaluated 
by using Maxwell equilibrium distribution for rrizo 

p{m,o) = (22) 

where k is the normalizing constant. The typical ex¬ 
pected value at room temperature is around rrizo ~ 0.99 
{Oq ~ 0.5°).!^ To compare the switching times it is im¬ 
portant to remember that conventional and SHE devices 
differ in two aspects. On the theoretical level, in conven¬ 
tional devices switching occurs at aj ^ a, while aj ^ 1 is 
required for SHE switching. On a practical level conven¬ 
tional devices can bear smaller currents due to heating 
problems. Thus achieving aj ^ 1 in them is problematic. 
In view of that, we perform two comparisons. First, we 
compare the SHE and conventional switching times for 
aj = 0.5 and h = 0. Here we get tp ~ 14 and ts ~ 5. 
Given the same normalized spin torque strength, con¬ 
ventional device is faster than the SHE one. Second, 
we compare the two devices operating at their critical 
switching current with a small field, say h = 0.02, point¬ 
ing toward the — 2 : direction. For the SHE device we 
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again use = 0.5 and the resulting switching time does 
not change much, tp 13.5. For a conventional devices 
we use aj = a, then tg ~ 29/(a. In this sense the SHE 
switching turns out to be much faster. In addition, since 
the initial condition is a statistical average, the switch¬ 
ing time estimated in this fashion may cause an non- 
negligible error in experiments. 


VI. SUMMARY 


The method of Refs. |29] and |30| provides a framework 
that can be applied to find the critical switching surfaces 
for any magnetic system with arbitrary anisotropy and 
spin torque in an exact fashion. In this article we calcu¬ 
lated the three-dimensional critical surface for a SHE bi¬ 
layer system with perpendicular anisotropy and in-plane 
current. Eor external fields in the x — y plane, the SHE 
induced spin torque not only shifts the existing equilibria, 
but also generates two new off-meridian equilibria. 

Firsts we derive an analytic formula for the stability 
boundary for the off-meridian equilibria — in previous 
numeric research this boundary was not distinguished 
from the boundary their existence. Second^ in contrast 
to the other authors we discuss the switching phase dia¬ 
grams in the field space at a constant current, and plot 
the modified astroid and the oval-shaped stability region 
of the off-meridian state. We then show how our qualita¬ 
tive description of the evolution of the constant current 
switching boundaries can explain the results of the other 
authors obtained for variable currents. Thirds we discuss 
in detail the evolution of equilibrium points and the char¬ 
acter of their destabilization on the switching boundaries. 
This allows us to put forward a qualitative understanding 
of the complicated hysteresis processes that are found in 
SHE devices. Fourth^ we point out that while a large cur¬ 
rent is required to set magnetization into the off-meridian 
state, it remains in this state when current is decreased 
to values that are a times smaller, and thus can be stud¬ 
ied at low currents. Fifths we show that in the presence 
of stable off-meridian state (current-dominated regime) 
switching between “up” and “down” equilibria happens 
through the off-meridian state. We consider an example 
of such two-stage switching and find the relationship of 
the corresponding switching time with that in conven¬ 
tional spin-torque device with magnetic polarizer. Here 
we find that, depending on the current limitations not 
related to spin torque physics, either of the two devices 
can operate faster. 
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Appendix A 


The explicit expression of each component in Eq. © 
can be derived as 

def = - dees + dghsp^, (Ala) 

def^ = -d0 + dehsp'^, (Alb) 


0 cos 6 » fir. , 

smu \smu 


1 


s ) , 


sin^ 
cosO 

{d0e - /isp®) 

sm6' ^ ^ 

-aLe + 


sm 6> sin0 ^ 


(Ale) 


(Aid) 


Appendix B 


Eor a planar linear system of the form X = AX, the 
eigenvalues of the 2x2 coefficient matrix A can be cal¬ 
culated in terms of its trace and determinant aPl 

/i-t = ^ ^trA ± y^(trAp'^^^detA^ . (Bl) 

Therefore knowing trA and detA tells us virtually every¬ 
thing about the geometry of its solutions. 

Besides stability, the types of an equilibrium is also of 
importance in understanding the switching process. An 
equilibrium of the same stability (except saddle point) 
can be a no de or a focus, depending on whether the eigen¬ 
values (Bl) are real or complex. Therefore the differen¬ 


tiation of the focus set and the node set requires another 
pair of critical values Ac± which satisfies 


(trA)^ - 4detA = 0. 


(B2) 


It can be demonstrated that the relationships Ac+ ^ A+ 
and Ac- ^ A_ are always satisfied. 

The classification of stability and equilibrium type in 
terms of the eigenvalues and A are summarized in TABLE 
|T]and in TABLE [TTj respectively. The complete dynamic 
analysis of an equilibrium needs to take into account the 
two factors together. 


Appendix C 


We mentioned the equilibrium types of the two 
off-meridian equilibrium states and their evolution in 
Sec. |IV A| and in Sec. |IVB 2[ To quantitatively under¬ 
stand the evolution of these two states, we also need to 
find the cr itica l A which separates nodes and foci, i.e., to 
solve Eq. (B2). It’s solution in the off-meridian case is 
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TABLE L Classification of stability 


Stability Type 

Equilibrium Type 

Eigenvalue Equivalent 

trA-detA Equivalent 

A Equivalent 

Sink 

Stable focus 

Stable node 

Complex Re[p,±] < 0 
Real fj,- < < 0 

trA < 0, detA > 0 

f A > Max(AT, A+) 

I At < a < A 

A W 

> > 

1 1 

Source 

Unstable focus 
Unstable node 

Complex Re[p,±] > 0 
Real 0 < fj,- < fj,^ 

trA > 0, detA > 0 

J A < Min(ATjA_) 

1 A+ < A < At 

if At ^ A-(- 
if At > A-(- 

Saddle 

Saddle point 

Real fj,- < 0 < p,-|- 

detA < 0 

A < A < A+ 


TABLE IL Classification of focus and node 


Equilibrium Type 

Eigenvalues 

trA-detA Equivalent 

A Equivalent 

Node 

Real 

(trA)^ — 4detA > 0, 

^etA > 0 

Ac- < A < A- or A+ < A < Ac+ 

Eocus 

Complex 

(trA)^ — 4detA < 0 

n 


A < Ac- or A > Ac+ 


^ Note that (trA)^ — 4detA < 0 guarantees detA > 0. 


given by 

Ac± — cos^ ^ ^ \ 

-^ajasinO coscj). (Cl) 

The two critical values Ac± give another two surfaces 
in the parameter space, one above A+ and the other be¬ 
low A_. Since both the two off-meridian equilibria have 
A = — cos^6> < A_, we only need Ac- to determine the 
equilibrium type. 

The two off-meridian equilibrium have different equi¬ 
librium types. For the one with 7r/2 < 0 < 37r/2, we 
have A_ ^ A ^ Ac-, therefore it is always an unstable 
node. The other one which satisfies —7r/2 < (p < 7t/2 


may change the equilibrium type as field changes, we 
can find two critical curves by equating A to A- and to 
Ac-, respectively. The former gives the destabilization 
boundary SqMi the latter corresponds to the type tran¬ 
sition boundary Sc of the analytic form 


(Vq;2 + 1-1) / ^ 

4a2a| ( ay' 

(C2) 

The transition boundary Sc touches Sqm but never 
crosses it. It separates the unstable region of the off- 
meridian equilibrium into node and focus regions, as 
shown in Fig. 
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regions large enough to be seen. 
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